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Abstract 

We analyze the solutions of the Schrodinger equation with the low frequency initial data 
and a time-dependent weakly random potential. We prove a homogenization result for the low 
frequency component of the wave field. We also show that the dynamics generates a non-trivial 
energy in the high frequencies, which do not homogenize - the high frequency component of the 
wave field remains random and the evolution of its energy is described by a kinetic equation. 
The transition from the homogenization of the low frequencies to the random limit of the high 
frequencies is illustrated by understanding the size of the small random fluctuations of the low 
frequency component. 


1 Introduction 


We consider the Schrodinger equation 

idt4>(t, x ) + - A (f>(t, x ) — eV ( t , x)(j>(t, x) = 0 (1-1) 

with a low frequency initial condition of the form 


<K 0,x) = 4> 0 {£ a x), (1.2) 

with some a > 0. Our goal is to analyze the long time behavior of cj)(t,x), and understand the 
energy transfer from the low to high frequencies that comes about from the inhomogeneities in the 
random media. 


We define the Fourier transform of / as 

f(0 = [ e~^' x f(x)dx, 

JR d 


and assume that V (t, x) is a stationary mean-zero Gaussian random held with a spectral represen¬ 
tation 



(27r) rf 


(1.3) 
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Here V ( t, dp) is the stochastic measure and V(t, dp) = V*(t, —dp), so V is real-valued. Its covariance 
function and power spectrum are 

R(t, x) = E{H(s, y)V(s + t, y + x)}, R(u),£)= [ R(t,x)e~ lut ~ l ^' x dtdx. 

The spatial power spectrum (the Fourier transform of R(t, x) in x only) has the form 

R(t,£) = f R(t,x)e~* x dx = e-°®WR(£), (1.4) 

Jm d 

where R(£) € L 1 (M d ), and the spectral gap 0(C) > 0, so that 


RM 


m)m 

w 2 + 0 2 (C)’ 


(1.5) 


By Bochner’s theorem, we have R,R> 0. Throughout the paper, we assume that 

€L 1 (R d )nL°°(R d ). 

0 (P) 


( 1 . 6 ) 


The compensated wave function 

The standard approach to an understanding of the behavior of the solutions of the weakly random 
Schrodinger equation is in the context of the kinetic limit OEa sea El ElE], through the study 
of the Wigner transform of the solution (the phase space resolved energy density) [9j. Our work 
here is closer in spirit to [H [TO] that focused not on the weak limit of the energy density of the 
solution but on the strong limit of the wave field itself. In order to motivate the “correct” way to 
this end, let us mention that after a long time the phase of the wave field acquires a large factor: 
for instance, setting 0 = 0 in (11,11) leads to an explicit expression 

0(t,C)=e-^/ 2 0( 0,0 

for the Fourier transform of the solution. Thus, a convenient object in the context of long time 
behaviors is the compensated wave function 

^(t,0=e <l€|2t/ %0, (1-7) 

which eliminates the deterministic component of the phase. This procedure is also known as phase 
conjugation in the engineering and physical literature. The surprising miracle is that after this 
simple-minded phase compensation, the wave field has a non-trivial limit. 


Loose end #1: the high frequency initial data 

We first describe the results of [4] obtained when the initial data for (11.11) is not slowly varying: 

^(0, x) = 0000), 
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that is, a = 0 in (11.21) . Let us set 


D(p,0 = 


2 R(p) 


(2^[ 0 ( P )-i(iei 2 -ie-p| 2 )/2] 

It is straightforward to check that 

(27r)V(p) + (iei 2 -|e- P| 2 )74] i^Y 

One of the results of [3] is that if 


, D(0= / D(p,Qdp. 


1 *(^^). 


( 1 . 8 ) 


(1.9) 


R(p) 

dip) 


£ L\R d ), 


then on the time scale t ~ e 2 , the compensated wave function corresponding to the initial data 
with a = 0 converges pointwise in distribution to a Gaussian random variable: 


- f i|£| 2 t - 1 , 

<%,£)^ =*MZ)e~* mt + Z(t,0- 

Here, Z(t , £) is a centered, complex valued Gaussian with the variance 

E{\Z(t,0\ 2 } = W(t,0 - \MH)\ 2 e- ReDi0t - 

The function W solves a (space-homogeneous) kinetic equation 

d t W= [ R( 

J«L d 

with the initial condition 




( 1 . 10 ) 


( 1 . 11 ) 


( 1 . 12 ) 


w(o,t) = \MH)\ 2 - 

This result is consistent with the aforementioned “traditional” kinetic equation approaches. 


Loose end #2: homogenization of the very low frequencies 


The results in the high frequency regime (a = 0) should be contrasted with the analysis of Bal and 
Zhang in mm for the case a = 1 in <tm performed for time-independent potentials. For the 
initial value problem 

i(j) t + ^ Acp - eV(x)cf> = 0, (1.13) 

4>(0,x) = (j)o(£x), 


with a mean-zero Gaussian random potential V(x), they have established a homogenization result: 



3 










converges in probability, as e —>• 0 to a deterministic limit cf>(t,x), which satishes the Schrodinger 
equation 

i(f>t + - Vcj) = 0 , ( 1 - 14 ) 

4>(0,x) = (f> 0 (x). 

The effective potential is constant and is given by 

- f R(p)dp 
.L \p\ 2 ' 

Let us mention that the choice a = 1 is special, as then the overall phase of the solution at the 
times t ~ e~ 2 is 

^e 2 |e| 2 = 0(l), 

so that no phase compensation is needed. 

Homogenization of the low frequencies 

Summarizing the above results, while solutions of ( 11 . 11 ) with the high frequency initial data have 
a random limit on the time scale t ~ e~ 2 , as in (11.10[) . solutions with the “very slowly varying” 
initial data as in ( 11 . 131 ) are homogenized on this time scale - their limit is deterministic. The first 
goal of this paper is to understand where the transition between the two regimes occurs - this is 
the motivation for introducing a general a > 0 in ( 11 . 21 ) . It will turn out that the homogenization 
result (formulated for the compensated wave function) holds for all a > 0 - that is, no matter 
how “relatively high” the low frequency of the initial condition is, solution has a deterministic 
limit at times t ~ e~ 2 . However, we will see that, unlike in the setting of [HE], the temporal 
fluctuations of the random potential lead to an effective potential with a non-trivial imaginary part. 
This means that the homogenized field loses mass in the limit. This loss of mass is attributed to 
the energy transfer to the high frequencies, which, as we show, account for the mass missing in the 
low frequencies, do not homogenize, and satisfy a kinetic type limit. We also analyze the random 
fluctuations of the low frequency component of the wave field and characterize the corrector to the 
homogenized limit. 

More precisely, we consider the Schrodinger equation 

idt4>(t, x ) + -A cf>(t, x) — eV ( t , x)cf)(t, x) = 0 (1-15) 

with a low frequency initial condition 


4>(0,x) = 4>o(nx), (1.16) 

with k<C 1. The Fourier transform of the initial condition is 

$( 0 ,£) = K~ d cf> 0 (-). 

Kj 
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Thus, if the function 0o(£) is of the Schwartz class, 0(0,£) is concentrated on the wave vectors £ of 
the size 0{k). While the Schrodinger equation with a time-dependent potential conserves the total 
mass: 

M(t)= f \fi(t,x)\ 2 dx = f |0(O, x)\ 2 dx, (1.17) 

J R d J R d 

the total energy 

E(t)= [ [|V0| 2 + eR|0| 2 ]oh (1.18) 

J R d 

is not conserved, unlike for time-independent potentials. Thus, even if the mass is initially concen¬ 
trated in the low wave numbers, after a long time evolution it may spread to 0(1) frequencies as 
well. As the potential is weak, the time it takes for the mass to spread over a range of frequencies 
will be long. 

We consider the long time behavior of the solution, on the time scale of the order t ~ e -2 , 
when the effect of the weak random potential will be non-trivial. We will first consider the “low 
frequency” rescaled compensated wave function: 


, ^ f. in? 

= K d $(- 2,f<)e 2^ (1.19) 

with the initial data 0 e (O,£) = 0o(0- This allows us to study the low frequency component of 
the solution - wave numbers of the order O(n). A straightforward computation shows that this 
function is a solution of the following integral equation 


, A , 1 f [ V ^’ a P) *K 2 (l€| 2 -k- £ | 2 )^r lit P\J 

^e{t,0=M0 + — / s d e “ (s,£ - -)ds. 

is J n J R d (2 tt) a K 


rt r v (-p’ d P) ki 2 -k- £ i 2 ) 

l£ Jo Ju d (2tt 

We have the following result for the low frequencies. 

Theorem 1.1. Assume that k = s a with a > 0. Then, for fixed t > 0 and f £ W l 

—> 0(t,£) = 0o(O e_ ^ D( '°' )f probability as s ^ 0. 


( 1 . 20 ) 


( 1 . 21 ) 


Let us stress that £ = 0(1) in the argument of the function ?b e (t, fi) corresponds to £ = O(k) in 
the argument of the function 0 - Theorem 11.11 addresses the evolution of the low frequencies of the 
solution of the Schrodinger equation with a slowly varying initial condition. Recall that 

0(0) = [ 7 — , 2R } pS> , , 9 . . dp, (1.22) 

Jr* (2tt) q! (0(p) + i\p\ 2 /2) P 

and, as g(p) > 0, we have ReO(O) > 0. Therefore, the passage to limit e — > 0 in (11.2111 induces a 
loss of the L 2 (M . d ) norm: while 

||0e(^) ') ||l 2 = II^oIIl 2 ) 

as can be seen simply from the definition of 0 e (t,£), we have 

ll^(i,-)IU 2 = ll4llL 2 e- ReD ^/ 2 <||0o||L 2 - 

The natural question is how does the loss of mass happen, and where does the mass go? Mathe¬ 
matically, there is no contradiction, as we will show the convergence in Theorem 1 1.1 1 is not uniform 
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with respect to £ £ M. d . From a physical point of view, as we have mentioned, the time dependence 
of the random potential breaks the conservation of the energy (11.181) . which allows the mass to 
escape to the high frequencies. Let us mention that in the time-independent case [2], where the 
conservation of the energy prevents the escape of mass from the low frequencies, it is shown that 
the mass is conserved as well. 


Generation of the high frequencies 


We now investigate the generation of the high frequencies in the above setting. Once again, we 
consider the solution <p(t, x ) of (11.151) with the initial data (|1.16l) . We stress that in all our results the 
initial condition (11.161) is the same - various rescalings in Theorem 11.11 above and Theorems 11.21 [L3l 
and 11.41 below correspond to zooming into various frequency ranges in the same solution. Our next 
goal is to understand how the mass escapes from the low frequencies (those of the initial condition) 
to the high frequencies, generated by the interaction with the random potential. As we are now 
interested in the high and not the low frequencies, we define the compensated wave function not 
quite as in (II. 19|) . but as 

rl ~ f 

(1.23) 


/ x d ~ , t . illlt 
= K' 2 4>(^2,C)e'^ T ~, 


so that the frequency is not rescaled. The initial condition for is 

®e(0,0 = K~ d/2 M£/K). 

The pre-factor K d / 2 in (11.231) is chosen so that we get a non-trivial limit. This function solves the 
integral equation 


^(*>0 “ Zdj2 + - 


V(jl, d P) J(|^|2-|g- p |2) 


J o jRd (2ir) d 


^®e(s,£ ~ p)ds. 


(1.24) 


The following result explains the loss of mass observed in Theorem 11.11 and tracks the generation 
of the high frequencies. 

Theorem 1.2. Assume that n = e a with a > 0, then for fixed t > 0 and £ ^ 0, we have 

^(LO =>- Z(t,£f) in law as e —> 0, 

where Z(t,f) is a centered, complex valued Gaussian random variable. Its variance Ws(t,£) is the 
solution of (11.121) with the initial condition 1^(0,£) = ll^oll!^)- 


The variance Ws(t,f) can be explicitly written as a series expansion 

w s (t,o = w s , b (t,z) + Ws,s{t,o, 

with the ballistic part 

Ws, b (t,Z) = ||^o|| 2 e- ReD( 0 )t « 5 ( 0 , 


(1.25) 
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and the scattering part 

OO n n k 

w s ,s(t,t) = Y UoWi / dv I <ir\[c <?* ... /',) 

k=1 J 0=v k+1 <Vk<...<vi<vo=t JR kd j_ g 

k 

x J] R eD(Pj,£ - ... - Pj-im ~ Pi - ••• - P k ). 

3 = 1 

Let us mention that Wfit,fi) = W^Jt. f) when £ / 0, that is, only the scattering part contributes 
to the variance in Theorem P We also observe 

[ W s , b (t,Z)d£= UoWle-^ 0 ^, 

JM. d 

which equals to the mass of the low frequency waves. 

Theorems o and P describe the dynamics of ED on different scales of the frequency do¬ 
main. In the former case, the low frequencies are zoomed in, and we find a deterministic evolution 
(homogenzation). In the latter, we track the high frequency component of the solution, so that the 
low frequency initial condition shrinks to a point source at the origin, which generates the high 
frequency waves. 


The fluctuation analysis in homogenization regime 


We now return to the analysis of the behavior of the low frequencies. According to Theorem 11,11 
the compensated wave function homogenizes for the low frequencies, hence the next interesting 
object is the fluctuation, which we define as 

Ue(t,£) = -ijWfefoO -E{^e(t,£)})- 

hi ' 

Here, ip E (t,£) is defined as in (11.191) . Heuristically, since the homogenization limit in Theorem 11.11 
captures the ballistic component of the wave held, we expect small random fluctuations consisting of 
the remaining scattering components. Indeed, we will see that the fluctuation exhibits a kinetic-like 
behavior. Let us set 


0 if a E (0,1), 

W a (t, £) = -D(0, 0)e~ D< '°' )t J <£o(f - p)M£ + p)e~ M ‘ 2v dpdv if a = 1, 

—D{ 0, 0)te _D ( 0 ^ f 4>o(£,- p)4> o(£ + p)dp if a > 1. 

, J R d 


(1.26) 


Theorem 1.3. Assume that k = e a , then for fixed t > 0 and £ E we have 


Ue(t,0 => Z s {t,£) = X s (t,g) + iYs(t,£) as e ->• 0, 
where X$,Y$ are centered, jointly Gaussian random variables such that 

E{\Z 5 (t,0\ 2 } = W s , s (t,0), 
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and 


E{Z 5 (i,£) 2 } = W a (t,£). 


Therefore, we can write 

4(U) = ER(U)} + «' /2 w £ (U), 

and Theorem P shows that when k = e Q , with a < 1, the fluctuation £Y e (f,£) is approximately 
distributed as Zg(t, 0), a centered complex Gaussian random variable with variance Wg yS {t, 0). This 
is similar to the result of Theorem P for the high frequency, albeit the variance is now given by 
the transport solution evaluated at the origin £ = 0, since we are now in the low frequency regime. 
If we let a —> 0 (which is the same as k —>• 1, so that the initial condition is less and less slowly 
varying), then, formally, ^> e (i, £) is distributed as 

MO^ mt + z s (t, o), 

which is consistent with (11.101) . That is, Theorem 11.31 also interpolates between the deterministic 
limit for the low frequencies and the random behavior of the high frequency component of the 
solution. 

The Wigner transform of the random fluctuation 

Besides the pointwise fluctuation for a fixed £ E R d , we also consider the fluctuation of '(/y (t, £) as 
a wave field. The tool we use is the Wigner transform for some (3 > 0: 

I W,(t,x,i)=( + (1.27) 

J R d 2 2 (2vr)“ 

Let Wg be the solution to the kinetic equation 

dtW + £ • V X W = f R( |P|2 ~ |g|2 ,p - £)(W(t, x,p) - W(t, x, £))^, (L28) 

J R d 2 (27T) a 

with the initial condition 

^5(0, x,£) = ||^o||l<K£)<5(x), 

and Wg t b,Wg >s be the ballistic and scattering component of Wg, respectively: 

W 5)b (t,x,0 = ||^ 0 ||l ( l(£)<5(x)e- ReD ( 0 ) t , 

and 

W 5 , s (t,x,0 = jr\\M 2 2 [ dv [ dP Y[e-^- v ^ ReD ^-- p i) 

fc=1 do=v k+1 <v k <...<vi<v 0 =t J R kd j_ q 

k k 

x n ReD(Pj, £ - ... - Pj_ i)J(£ - Pi - ... - P k )8(x - £t + ^ P jVj ). 







Theorem 1.4. Assume that k = e a , a £ (0,1) and a+/3 = 2, then for any test function p £ S(M 2d ) 
and t > 0, 


f 

Jr- 


W e (t , x , £)<p* (x, £,)dxdf 


Ws !S (t,x,0)<p*(x,£)dxd£ 


in probability as e —> 0. 


As Theorem o indicates that the ballistic component of transport solution gives the low 
frequency behavior, we conclude from Theorems 1 1.31 and II. II that the small random fluctuations are 
described by the scattering component of the solution of the kinetic equation. 

This paper is organized as follows. First, in Section [2] we present the Duhamel expansion and 
the corresponding diagrammatic expansions and the moment estimates that are needed for the 
proofs of all theorems. Section [3] contains the proof of Theorem 11.11 Theorem 11.21 is proved in 
Section U Finally, Theorems 11.31 and 11.41 are proved in Section [5) 

Acknowledgment. This work was supported by an AFOSR NSSEFF Fellowship and NSF 
grant DMS-1311903. 


2 The Duhamel expansion and the moment estimates 

Theorems 11.1111.2111.31 and 11.41 are all proved using the moment method. For the convergence 

V*e(</>o(£) e ^ D ^ t , 

in probability (Theorem II.ID . it suffices to show the convergence of E{if £ (t, £)} and K{\ip e (t, £)| 2 } 
to their respective limits. For the convergence in law of 1 F e (t,£) and U £ {t,(f) to a Gaussian in 
Theorems P and 11.31 respectively, we need to show the convergence of the corresponding mo¬ 
ments E{'F e (t, £)) W } and E{U s (t, tf) M (Uf{t, £)) N } for any M. N £ N to their respective 

limits, which makes the analysis slightly more computationally heavy. In this section, we perform 
the preliminary moment estimates that are needed in the proofs of the theorems. 


The Duhamel expansions 


All moment estimates rely on the Duhamel expansions that we now recall. From now on, we 
will set k = £ a . For the low frequencies, we can iterate the integral equation (jl.20p for the 
function if e (t,£), and write the solution as a series 


with the individual terms 

1 f 


fn,e{t ,0 = 77 


(iz) n JA n (t) 


n=0 




Pl + . . . + p n 

£ a 


), 


( 2 . 1 ) 


( 2 . 2 ) 
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and the phase factor 


(2.3) 


Gn(C,S (n) ,p (rl) ) = J^(|£-Pl - ••• -Pfc-l | 2 - |f -Pl - ••• ~Pk\ 2 )'^- 

k =1 

Here, we used the convention /o, e (t,£) = $o(0> an d have set po = 0, = (pi,... ,p n ), as well 

as s( n ) = (si,..., s n ). We have also defined the time simplex 

A n (t) = {0 < S n < . . . < Si < t}. 


For the high frequencies, the solution \I/ e (t,£) to (| 1.241) is similarly written as 


with 




n=0 




T otd/2 


M 


£ ~ Pl ~ ■ ■ ■ ~ Pn ■ 


(2.4) 


(2.5) 


and 

^°’ £ ~ e ad/2^°^^a)' 

The key “bureaucratic” difference between the Duhamel expansions (12.21) and (12.51) for the func¬ 
tions ip e (t,£) and 'F e (t,^) is that H > £. This will make the limits very different. 

The following lemma ensures that the solutions given by 40 and (12.4p are well-defined and 
we can interchange the summation and the expectation when computing the moments. Its proof is 
exactly as that of [H Proposition 3.8]. 

Lemma 2.1. Fix £ > 0, M, N € N. Let g n , £ = / n , e or F U}£ , then 

l®{5mi,e ■ ■ ■ 9m M ,e9ni,e • • • 9n N ,e }| < C £ (mi,.. .,m M ,ni,... ,n N ) (2.6) 


with 


OO 

E 


Y C £ (mi,... ,m M ,ni,... ,n N ) < oo. 


mi,...,mM=0 ni,...,rijv=0 


The pairings 

Now. we discuss in detail the calculation of the moments 


where g n , £ = f H)£ or F n>e , and 


• • • 9m M ,e9n 1) £ ■ • • 9n N ,e}i 
M N 

Y / rn i + 'Yrij = 2k, 
i= 1 j =1 
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for some k € N (if the sum is odd, then the moment is zero by the Gaussian property). We have 


^{9mi, S ---9m M ,e9n 1 ,s---9n N , E } = 0 £ ) E ^ =1 Uj 

x / dsdu / E{/M,Af} 

JAm, (t)x...xA„„(t) 




A/ TV /o 7) 

h,M,i /ijVj J 

j=l 


2 — 1 


with 


Im,n 


i2 kd 


n^,d P1>1 )... vc-^,d P1 , mi )... v(^,d PMtl ) ...vc- 


-,dpM, 


(2?r 

x V*(^M jdgi>1 )... ,dq 1>ni )... V*(^,dq Ntl )... V*(^^-,dq N>nN ), 


m M ) 


and the phases 


M 


5m = £ G m („, •f'GI’" 0 )/* 2 , 5»=E I.'”*’)/® 2 , 


N 


(ru) (rij) \ 


2—1 


2—1 


with 77 = or £, depending on whether g ne = / n e or F ne . The initial conditions appear as 

u 2 it Pi,l + • • • + Vi,mi N u * 2*It + • • • + Qj,nj x 

h M ,i = Ms --)> h N,j = Ms --) 


when g n , e = f n>e , and as 
hM,i = £~ ad/2 M 


£ Pi,l ■ ■ ■ Pi,mi 




i 1j ■ ■ ■ Qj,r. 


when g n , £ = F nfi . 

Using the rules of computing the 2k— th joint moment of mean zero Gaussian random variables, 
we obtain 

E {Im,n} = E II (27r)- d e-®^)l’"-^l/ e2 <y(«; l + w r )R{ Wl )d Wl dvo r . (2.8) 

F (vi,v r)£F 

The summation is taken over all allocations of the set of the vertices 


{^1,1 j • • • i 'Sljmi) ■ ■ ■ j SM,1, • • • j $M,mM ^1,1 j • • • > ^l,ni j ■ ■ • j UN,li • • • i ^N,iin }■ 

into k (unordered) pairs (recall that J2iLi m i + Ylf=i n j = 2fc). We call each allocation a pairing. 
In f|2.811 . vi,v r are the two vertices of a given pair, and wi,w r are the respective p, q variables, that 
is, wi = pij if vi = Sjj and wi = —qi,j if vi = Uij. The same holds for w r . We will also write a pair 
as an edge e = (vi,v r ). Note that the order of vi, v r does not matter here since g, R are both even. 


A uniform bound on the pairings 

We recall the following general bound. 
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Lemma 2.2. Let g n>£ = f n>e or F n>e , then we have, for all e £ (0,1] 


|E{<7mi,e ■ • • 9m M ,e9ni,£ • • • 9n N ,e \I — 


(2k — 1 )!! 


a=iK)!nf=iK)! 


c k 


(2.9) 


with some constant C depending on t, £, i?, g. Tfere 2k = X)f=i + ^27=1 n j- 


iJV 


Proof. The proof is close to the case g = f n>£ and a = 0 which is already contained in [3]. 
We present it, together with the required modifications, for the convenience of the reader. By 
symmetry, the RHS of (12.711 can be bounded by 


M N 


TLm~, — Jt—tv , W 7 T / dsdu [ |E{ W |fll^lIll^|. ( 2 -! 0 ) 

nzi(mi)\nUM- £ jR2kd 11 7=i 


In the case when g n)E = / n)£ , we bound 

M N 


iNlIft 


< 


M+Af 


jl — liv'uiioo > 


*=i i=i 

then for a given pairing T, we have 


-L- [ dsdu [ IT (27r) d e Vr ^ £2 5(wi + w r )R(wi)dwidw r < C k , 

e J[0,t] 2k JM?kd ■ L ^ 

L ’ J (VI,Vr) G-F 

where we used the integrability of R(p)/q(p)- Thus, ( 12 . 101 ) can be bounded by 

yofcn? II M+N _ (2k— 1)!! k ~ , { m+N 


rrM ( \\y\N / N|^ M 


oM / unN / ll^ollc 


( 2 . 11 ) 


1 , 7=1 \JJ j- 

In the case when g n ^ = F n>e , we consider £ 7 ^ 0 and integrate w r and bound (12.101) by 


1 




yyr / * <f “ / E n 


M N 


,-Q{wi)\vi-Vr\/e‘ 


R(wi) \hM,i\ n 1^1 n 


duii 


d ■ 


I .7 — 1 V ,(/ 7 . 

J [C ),t] 2k R kd 

For a given pairing F, we have 


F (yi,v r )eT 


4 4 (2ir. 

i =1 j=l Z=1 v ' 

( 2 . 12 ) 


\h M ,i\=e ad/2 \M{^;), \h* N J=e Qd/ 2 |^o|(^-), 


(2.13) 


where 


Fi f, Pi, 1 • • • Pi,m,ii Qj f, Qj,l ■ ■ ■ 9j,np 

subject to the conditions 

wi + w r = 0 when (vi,v r ) £ T. 


(2.14) 
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The difference with the previous case are the factors £~ ad / 2 in (12.1311 . Note that if Pi = £ or Q j = / 
(this may happen because of (12.141) 1. as £ ^ 0 is fixed and <j>Q is rapidly decaying, we may simply 
use the bound 

e~ ad/2 \M {^)<C- 

For i, j such that Pi, Qj ^ £, to deal with the large factors in (12.131) . we change variables as follows. 
Take some % with Pi ^ £, so that 

Pi, 1 T • • • T Pi,mi 0. 

We pick any variable p from {pi,i, ■ ■ ■ ,Pi, mi } (note the number of elements here can be strictly 
smaller than m* since we have already integrated out the variables w r ), and change p to p' = Pi/e a . 
The variable p = wi was paired to some pj or q.j = w r as in (J2.14D . Thus, after the integration of 
w r , p' will also appear in a unique fiM,i which equals to some Hmj or h* N j ■ We use the bound 

\h M ,i\ <e~ ad/ 2 C. 

Thus, after the change of variable and taking into account the Jacobian of the change of variables, 
we have, with a slight abuse of notation 

I h M ,ih M ,i\dp < e~ ad/ -\(j) 0 \{p')e~ ad/ 2 Ce ad dp' = C\(j) 0 \{p')dp'. (2.15) 

Since the change of variable only relates to pi, all other Hms, h* N ^ are not affected. We continue 
the procedure, integrating out the p-variables one by one. If we are left with a single 

\h\ = e~ ad / 2 \(f> 0 \(Pi/£ a ) or £~ ad/2 \4>o\(Qi/£ a ) 

in the end, we change variable similarly, and estimate this term, together with the Jacobian as 

£ ad/ ~\(t>o\{p')dp < | 0 o| (p')dp'. (2.16) 


Overall, this change of variables will involve M + N momenta, and will eliminate all factors hM,i 
and h* N ■, and we will be left with an expression of the form 


/ dsdu [ TT ( 2 

J {vi,Vr)^ 

/ dsdu / TT (2 

J\o,t] 2k JR kd , 72^ 


M 


N 


'{0,t] 2k 

(jM+N 


n)- d e -^ wi) \ vi -^\/ £2 R{ Wl) n i h M ,i\ n \h* N jd 

(yi,v r )£F i= 1 j =1 


W 


£ 2 k 


TT )~ d e -^ Wl )\ Vl - Vr \l £2 


R(wi) JJ (27r)- d e- 9 ^ v ‘- Vrl/£2 R(z l ) 

(vi,v r )€iF :j 


(vuv^GiFi 


x II \M (wi)dw. 

(vi,V r )eT2 

(2.17) 

Here, (vi,v r ) G P\ denotes the pairings in which the momenta do not participate in the change 
of variables and ( vi,v r ) £ Ti denotes the affected pairings. The explicit form of zi that appears 
above is not important, so we do not specify them. The bounds (|2.15l) and (|2.16l) mean that the 
“participating” wi give us the factor 


II l^olM 

(vi,V r )&iF2 
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that appears in the last line of (12.171) . 

Next, we integrate in time. This brings about the product 

^ k TT R( w l) TT R( z l) 

( ft PF bm , ^f PF 0(*i)' 

\yi,v r )€-T\ \yi,v r )€.T2 

Using the fact that R(wi)/ q(wi) is integrable for the vertices in J-\. and that R(zi) / g{zi) is uniformly 
bounded for the vertices in 7^, we may integrate out all the momenta variables, showing that ( 12 . 121 ) 
is bounded by 

m rik_ (2fe - 1)!! ck 


ufii(my.uU^y c a=iK)!nf=iK-)! 


(2.18) 


This finishes the proof. □ 

Lemma 12.21 ensures we can interchange the limit e —> 0 and the summation, since 


E E 


(2k — 1 )!! 


= 0 B1 ,.^=o Uiil K)! njliK)! 


N 




(M + N) 2k C k 

-C = 2_ -—- < oo. 


k =0 


2 k k\ 


An estimate on non-simple pairings 

Now we need to consider more carefully the contribution from different types of pairings. First we 
can decompose the temporal domain A mi (t) x ... x A njv (t) according to all possible permutations 
of {si,!,... u NtnN } and write 


■ ■ ■ 9m M t s9ni ,e • • • 9n N ,e\ 


M N 


E -Tuv- [ dsdu [ E {lM,NV GM e- lGN ft h W II h kv 

~ (ie)Ei=i m i Ja 2k (t) Jm. 2kd 7i TA 


(2.19) 


i =1 j =1 


where cr 2 fc(t) = {0 < u 2 fc < ... < v i < t} and o = {ui,... , u 2 fc} denotes all possible permutations 
of {si ) i,...Ujv, rajv } such that a 2 k(t) / 0. By (|2l8|) . 

( 2 . 20 ) 


WmM = E II (2^)- d e-s^)l^-^l/ £2 «5(^ + uv)£(«;,)du>idttf r , 

? {vi,V r )&T 

where J 7 are pairings obtained by computing joint moments of Gaussian. We can write 

' ‘ ' 9mM,s9n\,e • • • 9nj^,e] ^ ^ ^ ^ C) 5) 

o' J 7 


( 2 . 21 ) 


with 



1 

v—'iV 

(-ie)^i= inj 

sM\vi-Vr\/e 2 S ^ Wi + Wr )R( Wl )d Wl dw r e iGM 


e -iQN 


M N 


h-M,i H j 
*=i i=i 

( 2 . 22 ) 
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and the symbol g = / or F indicates the dependence of J^ li n *^ on g n>£ = f n>£ or F n>e . 

Given a permutation a, we say that T a is a simple pairing if V 2 i~i,V 2 i form a pair for every 
index i = 1,..., k. The next lemma shows that the overall contribution of the non-simple pairings 
vanishes in the limit e —> 0. 

Lemma 2.3. Let g n£ = f ne or F ne , then we have 

EE ((*,^,€,9) 0, as e -> 0. 

* T^T a 


Proof. When g Uy£ = f n . £ , this is proved in |U Lemma 3.6]. The proof for g n>e = F n>£ is 
similar, using the same change of variables as in the proof of Lemma 12.21 We do not provide all 
details - but just mention the main simple point: for T ^ F a . i.e., the non-simple pairings, the 
non-consecutive times are paired, then the overall contribution from the time integration of the 
exponentials e~^ Wl ^ Vl ~ Vr ^ £ is too small since \v\ — v r \ is too large. We write 

\ J m 1 ,...,n* N ( a i J T ,£,,F)\ 


< 


1 ,. M N 

~ 2 k dsdu / IT (2TTy d e~ s( - Wl ^ Vl ~ Vr ^ £2 5(wi + w r )R(wi)dwidw r Y\\h M ,i\T\\h* N 

£ j9 - 2kd (v lVr )eT 11 7=1 


and by the proof of (12.171) , we have 


1 4,..,n^. 

/ ~<m+n r r 

/ dsdu / TT (2v T)- d e-^ Wl) \ Vl ~ Vr \/ e R{wi) TT (2tt)- R(-) 

£ Jw. kd , ■ L r_ x . , ■ L r =7r 

' [VuVrjeT 1 ( Vi,V r )e7F2 

x II \7>Q\{wi)dw. 

(vi,Vr)eT 2 


Then, using the fact that R(p)/q(p) is integrable and uniformly bounded, we only need to follow 
the proof of [H Lemma 3.6] using the aforementioned observation that the time integration will 
bring about too high power of e because of the exponential in time factors. □ 


The vanishing of the crossing pairings 
By Lemma 12.31 we have 

lirnE {g mi , £ ...gm M ,e9*n 1 ,e ■ • ■ 9n N ,e} = So ^ ff) ■ 

cr 

Let us define sets 

A-i = {Si,i, • • • ,Si,m J, Bj = {Uj ! 1 ,.. .,Uj !nj } with i = 1, ..., M, j = 1,..., N. 
Given a pairing J>, we say 


(2.23) 
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interact with each other if there is an edge (vi,v r ) G T a such that vi G S\,v r G £ 2 , and we 
write Si G4 S 2 . We say they are connected if there exist other sets such that Si -H- ... -H- 
S 2 . Thus, for a given permutation < 7 , we may decompose {Ai,Bj : i = 1 

into connected components. For example, if all variables in ^4i pair inside Ai, then A\ itself 
is a connected component. If all variables in A\ and A 2 either pair inside the corresponding 
set or pair with variables in the other set, and we have at least one edge joining A\ and A 2 , 
then {A\,A 2 } is a connected component, and so on. We let N c (F a ) be the size of largest connected 
component corresponding to F a . The following lemma shows the permutations with more than 
triple interactions do not contribute in the limit. This leads to a Gaussian limit in Theorems 11.21 
and 11.31 


Lemma 2.4. We have 


and 


Y = 0. 

a\N c {F a )>P> 


Proof. We first consider the case g n>£ = f n ,e- For a given permutation < 7 , if N c (a ) > 2, we can 
find the sets 

Si,S 2 G {Ai,Bj : i = 1,..., M,j = 

such that 5i o £ 2 - Let e be an edge joining S± and £ 2 , and hs 1 ,hs 2 be the initial conditions 
corresponding to £ 1 , £ 2 , then we have 


s~iM-\-N—2 /• /• 

<- 77 - / dsdu / TT {2'K)~ d e~^ Wl ^ Vl ~ Vr ^ e 5(wi + w r )R(wi)dwidw r \hs 1 hs 2 \, 

£ U2kd (nirter. 

(2.24) 

where other factors of the initial condition (fio are bounded by Recall that when g n ^ e = 

f n ,e, we have 


hM,i = 4>o{£, — 


Pi, 1 + • • • + Pi,; 


h* N>j = m- 


Qj, 1 + ■ ■ ■ + qjj 


We can assume 


I foil = I<£o|(£- p) and |fo 2 l = l0o|(C — 

for some Pi, P 2 after integrating out w r in (|2.24l) . It is clear that Pi, P 2 7 ^ 0 since they both contain 
the wi variable corresponding to the edge e. Now we integrate w r and the time variables to obtain 

J id n ^yl^olte - h)|*,|({ - p)|dw, ^ 0 , (2.25) 

Rkd («|,S r )e^r 1 


as £ —> 0 by dominated convergence theorem. 
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Next we consider the case g U}£ = F Uj£ . The following estimate holds 


< 


-t r r M N 

2k dsdu JT ( 27 r)~ d e^^ Wl ^ Vl ~ Vr ^ £ 2 S(wi + w r )R(wi)dwidw r Y\\h M ,i\[\\h* N j\. 

Ja 2k (t) Jm^d ^ i=i , =1 


Recall that now 


h M ,i = e~“ d/2 M S Hl e „" Km ), h>„j = s-^U 
If N c (a) > 3, we can find 


£ Qj, 1 • ■ • Qj,r 


Si,S 2 ,S 3 € {Ai,Bj : i = 1,. .., M,j = 1,..., N} 

such that Si ++ S 2 S 3 . We pick two edges linking Si to S 2 and S 2 to S 3 , and denote them by ei >2 
and e 2 i 3 , respectively. We also denote the variables corresponding to 61 , 2 , 62,3 by wi >2 , m 2 , 3 - Let hs t 
be the initial condition corresponding to Si, i = 1,2, 3, then we have 

\hsi \ = g- Qd/ 2 |^o|( g £C Pl ) for some P h i = 1,2,3. 

After integrating out the w r variables, it is clear that Pi contains the variable w 12 , P 2 contains the 
variables wi 2 ,w 2 3 and P 3 contains the variable 162 , 3 . We do a similar change of variable as in the 
proof of Lemma [2721 First, we change wi j2 so that (£ — Pi)/e a Pi. Second, we change 162,3 so 
that — P 3 )/e a ^ P 3 . Then we have 

\h Sl h s ,hs,\ =e- 3ad/2 \M^^)M^^m^^-)\dw v2 dw 2 , 3 

with some z that does not matter to us, as we simply bound (fio(z) by C. Now, we only need to 
carry out the same change of variable as in the proof of Lemma 12.21 for the remaining h. In the 
end, we obtain 

(~iM+N F ad /2 r r 

\J e mi ,...,n h ^^,^F)\ < -- / dsdu / n (2n)- d e-°^-^ R( Wl ) 

x JJ (2vr)- d e-^ Vl -^ 2 R(z l ) J] \fo\(w t )dw, 

{vi,Vr)€.T a ,2 (yi,v r )e.J 7 lTl 2 

(2.27) 

where (vi,v r ) € F a ,i denotes the pairings which are not affected by the change of variables, 
and (vi,v r ) G T a , 2 denotes the affected pairings, and, as in the analysis of (12.171) . the precise 
expression for z\ is not important to us. Clearly, the RHS of (12.271) goes to zero as e —> 0 because 
of the extra factor £ ad / 2 compared to (12.171) . □ 
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Pairings for the correctors 

We now describe analogous estimates that are needed in the analysis of the corrector 


n=0 

with 


^n,e(*,0 = fn,e(t,£) - IE {f n ,e{t,Q} 

Jm. nd e £ 


{ie) r 


>A„(t) 


Sl Sn ' -7_iG„(e“$,s<"V">)/e 2 i ft _ Pi _ _ —)ds 

g-a ‘ ‘ ‘ gat ' ’ 


where 


V(si, • •• ,s n ,dpi,.. . ,dp n ) = (27r) ncZ JJ R(sj, dpj) - E{JJ dp,)} 

\j=l 3 =1 

Let us discuss in detail the calculation of moments 

£ -ad(M+JV)/ 2 E{ ^ i £ . jr mM)6 jr* i(e ... for m 4 ,,njv€N, 

with 


(2.28) 


M iV 

rij = 2 k. 

i= 1 7=1 

Similar to (12.71) . we can write 

-ad(M+7V)/2|j 7 r «-* \ 

t ■ ■ ■ ^m M ,e^ ni! £ ■ ■ ■ ^n N ,ei 


1 1 f 

Li=l m i (_i F )I2j=l n j .]/\ 


M N 


( is)T,i=i m i (_j e )Sj=i n i JA mi (t)x...xA„ N (t) 


dsdu [ E{1 m,nV Gm e~ iGN ft h M ,i if h *N 
j ^ kd t=1 J=i ’ 


where 


T ^l, 77 ll i 7 \ 7 . / .l 7 !, 1 7 7 \ 

AM,IV = L(^ 2 “: ■ • • ) g2 7 QPl.l 7 • • • j dPl,m\ ) ■ ■ ■ . . . , -^- > d PM,l, ■ ■ ■ ,dPM,m M ) 

_ -.,* 7 ^ 1,1 1 * 1,711 7 7 \ U N , n N , 7 \ 

X V ( ^2 ? • • • 5 ^2 ‘ ^5 • • • 5 ) • • • V ( ^2 j • • • 5 ^2 ‘ ? * * * i ^ QN,tin ) 5 


(2.29) 


(2.30) 


and 


hM,i = £ ad / 2 (i> o(£ — 


1 + ... +Pj- 


-), h* NJ =e~ ad / 2 m- 


Qj, ! + ••■ + Qj,r, 


e“ 7j £ u 

Previously, we have dealt with the expectation of a product of centered Gaussians. For Im,n, 
however, each factor V, defined in (I2.28P is a centered product of Gaussians rather than a product 
of centered Gaussians. The rules for evaluating the expectation of such objects are recalled in 
Lemma lA.il in the Appendix. Recall that we have defined the sets 


Ai = {s;,!,.. .,s itTn J, Bj = {u jt i,.. ■ ,u j>nj }, with i = 1,..., M, j = 1,..., N. 
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Given a pairing J 7 , we decompose {Aj, Bj : i = 1,, M, j = 1,..., N} into connected components 
according to the interaction between the s, u variables. Let N S (P) be the size of smallest connected 
component, then by Lemma lA.il we have 

E {1 m ,n}= Y n i^)~ d e~ s{wi)lvi ~ Vrl/£2 5{wi + w r )R{wi)dwidw r . (2.31) 

T:N a {T)>2 (vi,v r )eJ r 


In particular, it is clear that E {Pm,n} < ^{Im,n} and 

E {Im,n} - mM,N} = Y n (^)- d e- 0{wi)lvi - Vrl/£2 6(w l + w r )R( Wl )d Wl dw r . 

T:N a (T)=l (v^vr)^ 


Comparing (12.71) and (12.81) . to (I2.29P and (|2.31l) . we see that 

£ -ad(M+N)/2 E{ J? mi>£ . . . £ 

has exactly the same form as 



(2.32) 


• • • Fm M ,eF n i,e ■ ■ ■ ^njv.el (2.33) 

if we replace £ —> e a £ and impose the constraint N S (P) > 2 in (I2.33|) . Therefore, we can follow the 
same proof for Lemmas 12.2112.31 and obtain 


lim e -«*(W+JV)/2 E {^- 
€—^0 


mi,e ■ ■ ■ ^mMfi^ni.e 


n N 


} = 


V 

e —>0 

a:N s (T a )>2 


n m 4, 


r (a,J- ff ,e%F), 


where we recall J^ n ^ is defined in (12.221) . 

We should note that in the proof of Lemma 12721 for g n , e = F nj£ , we used the fact that £ ^ 0 so 
that 

e- ad/2 \Mi/n\ < C, 


and actually goes to zero as e —> 0. At this step, the analysis for & n>e can not proceed this way, 
as we have replaced £ i->- e“£. Instead, we use the condition N S (P) > 2, which implies that after 
computing moments, all the h factors in (I2.29P take the form 


h M ,i = e ad/2 0o{£ ~ ^), and h* N j = e ad/2 4> 0 {t, ~ ^), 

for some P,Q^0. If P or Q were to be zero, then Aj or Bj is not connected with any other set, 
which would imply N S {P) = 1. As P and Q are not zero, we only need to perform the same change 
of variables as in the proof of Lemma 12.21 

We may now follow the same proof as for Lemma 12.41 to obtain 


hm£- ad ( M+JV )/ 2 E {^ mi)£ 


• • ^rriM ,e 



Since N s (P a ) > 2 and N c (P a ) < 2, we have 


Y Yl J mi,.,n* N (^^F a C,F)- 

o:N s (F a )>2,N c (F a )<2 

(2.34) 


Ns(Pa) = N c (P a ) = 2, 


that is, all connected components corresponding to J- a contain two sets, which implies M + IV is 
even. 
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3 Homogenization of the low frequencies 

We now prove Theorem ll.il To show that 

ip E (t,£) —> in probability, 

we only need to verify the following result. 

Proposition 3.1. As e —» 0 ; we have 


and 


E{\Mt,0\ 2 } ^ \MO\ 2 e~ ReD{0)t . 

Proof. By Lemma 12.11 we have 

OO 

E (VVVVE^VV)}- 


n=0 


Lemma 12.21 ensures that we only need to compute 


limE {f n ,e(t,Q}, 
£—^0 


when n = 2k for some k G N. By Lemma 12.31 we have 

limE {f n ,e(t, 0 } = V hm J^(a, 

F —VO Z ^ P—Vf) 


£—^0 

It is straightforward to see that 

V<bV£,/) 

MO 

I us ! 

/— AiA lm>kd 

(vi,V r)&Fa 


V) 2fc Mm 


[ ds [ n (v 

J<T2k(t) JUL kd ’Tt-r- 


-d e -fl(to,)|«i-Vr|/6 2 _^^^ e i(|e“el 2 -|e a e-i«ll 2 )^|^ L 


2 e' 2 dw 


^o(0( _1 ) A 


r(p) d P \ 


W. IV (fl(p) + iH 2 )(2vr 


fc! 


and thus 


limE{^ e (i,f)} = V IimE{/ 2 fc >e (i,0} = E 1™ WV 0 f) = M0 e ^ (0)t , 

£ - * -* /=•-VII ' -* F -VII 


fc =0 


' £—^0 


k =0 


£—^0 


which is (13.111 . 
Since 


E{|W)| 2 } = E E {V(U)/^(U)}, 

m,n=0 


(3.1) 

(3.2) 


(3.3) 


(3.4) 
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by a similar discussion as in the proof of (13.11) . we have 


limE{/ m , e (t,0/n, e (*>0} = •?>>&/) 

£—^0 5 Z ' r— ^ 1 ) ’ 


£—>-0 


(3.5) 


In addition, Lemma 12.41 shows that 


Y J ™,n * (<T f ) = °> 

cr:AT c (cr)>2 


(3.6) 


so we are left with 

X] lin l J m,n*fo-7v /). 

(r:iV c (cr)=l 

However, N c (cr) = 1 implies there is no interaction between / m , e (f,£) and /^ £ (i,(), so m = 2k\, 
and n = 2/^2 are both even. The number of possible permutations is 

(Ah + k 2 )\ 

h \k 2 \ ’ 

and by the same calculation for (13.41) . we have 


Y J£ m,n* ( a - ^ f) ~ (li + k% 


11 „ fki+k 2 

-\MC)W-i) kl+k2 - 


cr.N c (cr)=l 


Therefore, we have 


=\M0\ 


h\k 2 \ l ' ruv " yi v *' (h + k 2 )\ 

2 (—f.D(0)/2) fcl (—£D*(0)/2) fc2 


(T)“ ( 


D 




k2 


2 ) 


Ah! 




(3.7) 


limE{|*(t,{)| 2 }= D ko(0| 2 < tP (°j /2)t ‘ ( iP 2” )/2) ‘ 2 - l^(0| 2 e~ RcD10) ‘, (3.8) 


ki,k 2 =0 


ki 


which is ()3.2I) . □ 


4 The high frequencies 

In this section, we prove Theorem 11.21 
Convergence of the mean 

We first show the convergence of E{T e (£,£)} for fixed t > 0 and £ / 0. 

Lemma 4.1. We have 

E{\h e (£, £)} —y 0 as e —>• 0. 
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Proof. By Lemmas 12.2112.31 we only need to show that 


lim J®(<7, = 0, 

£—>-0 


when n = 2k. It is straightforward to see that 

r n (a,T a ,d,F)=e- ad / 2 Mh 


(*) 2k 

Since £ ^ 0, we have 
thus 


l£) Jau(t) Jw kd , 


( Vl,V r )eJ r a 

£ -ad/20 o (£/ g a) 0 as £ —>■ 0, 

|J^(<T,7 r (T ,^F)| < Ce~ ad/2 4> 0 (£/e) ->• 0 as e ->• 0, 


and the proof is complete. □ 


Convergence of the variance 


Next, we look at the second moment. 


Lemma 4.2. We have 

E{|* £ (t,0| 2 }-> as e —>■ 0. 


The proof of Lemma 14.21 is very similar to [3, Proposition 3.12], and since Lemmas 15.11 and 15.21 
below follow the same blueprint, we will provide the details here for the convenience of the reader. 

Proof. By Lemmas I2.HI2.41 we only need to consider Jv,£, F) for fixed m,n £ N (in 

the present case, we automatically have N c {F a ) < 2). We write 

A — -(si,..., s m }, B — {ui , ■ ■ ■, 

with m + n = 2k for some k € N. According to the pairing J-„. {A, B} is decomposed into connected 
components. If A, B are “separate”, we have two factors of £ _ " rf//2 |0o|(C/ £a ) coming from the initial 
conditions, so by the same argument as in the proof of Lemma 14.11 we have 


n* (*L o j s ■ -F) —y 0 as s —y 0. 

Therefore, we only need to consider a such that A «->■ B. 

For a permutation a of A U B, the simple diagram J- a corresponds to 


AU B = f ...,v£,v k }, 


with 


A > v i 


>v£>v k 


and (yf,v i ) forming a pair, i = 1,..., k. Since A -H- B, there exists at least one pair such that vf 
and v^ come from different sets, and we call such pair a crossing edge between A and B. Assuming 
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the total number of crossing edges is N cr {a) > 1, the interval [0, t] is decomposed into N cr + 1 parts 
according to the position of those crossing edges, which we denote by 


’ i > rf > 


> r N > 
lvcr 


N c , 


with rf = vf for some j. and with the convention where r Q = t and = 0. We further 

denote by £ s j, i = 0,. .. ,N cr , the set of edges between the vertices r~ and rf +l that are of the 
form (sj,Sj+ 1 ), and by £ u j the set of edges between rf and rf, , that are of the form (uj,u,j+ 1 ). 
The corresponding sets of indices are denoted by 

A = {j ■ ( v f, v f ) € £ s,i,j = 1, • • • ,k}, 

and 

Si = {j ■ (v+,vj ) G £ Uji , j = 1,..., k}, 

with i = 0,... ,N cr . For a non-crossing edge (vf,vj), we denote Tj = 1 if vf,vj are s —variables 
and Tj — — 1 if they are u— variables. 

Recall that 


Jm,n* A At A A) 

= 77^7 . , n [ dsdu [ TT (2tt )~ d e~^ Wl ^ Vl ^ /£2 5(wi + w r )R{wi)dwidw r 

x g iGm (£,s (m) ,p (m) )/Pg -iG n (£,« (n) ,g (n) )/e 2 £ ~ ad l ( ^ ~ ^ >1 ~ ~ ( j Z ^ Z I ' ' ~ 

7UV /tOV ^.q, /’ 

where vi,v r are the vertices of a given pair, and wi,w r are the corresponding p,q variables, that 
is, wi = pi if vi = Si and wi = —% if ij = Uj. For a crossing edge (vi,v r ) = {rf , rf ) , the relevant p, g 
variables equal to each other due to 5{p — q), and we denote the corresponding wi = Pi, with the 
convention that Pq = 0- We also define Sj = 1 if rf is s —variable and s t = — 1 if rf is u— variable. 
With the above notation, we have 


J m,n* t,£,F) 

1 [ 

= —r—-— / dsdu 

M m {-ie) n J a2k{t) 



R(wj) 

(2ir) d 



Ncr 

pAiw 

3=1 


■~Pi- 



N C r 

n n 

j=0leAjUBj 


iTi(\£-...- p j\ 2 -\£-...-Pj- W i\ 2 )- 




^3wl^o|( 


i-Pi-...-p Nc . 


(4.2) 

Here, we have integrated out the variables w r in (14.11) . and changed the notation wi e-)- wj . To get 
rid of the extra factor e ~ ad , we change variables as before. Replacing 


Pn ct >-»■£ — -Pl — ... — Pn c ,—l - £ a PN cr , 
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and rewriting the terms in (14.21) associated with Pn ct using the new variable, we obtain 


Jm,n*( a iF<n£i F ) = T 




/ dsdu / dw 

Jao k (t) JRkd . ^ 


j:Wj^P Nc , 


R(Wj) - B ( W j) 3 j, 3 

{2ir) d 


- p i - ■ ■ ■ - P N cr -1 - ^ P N cr ) „-s {i-Pi---P Nrr -i-e a P Nrr .) 


r N cr r N c r 


(2vr) 


JVcr-1 

- n 

j =i \ 

Nci —1 / 

* n n 

j =o 

x |0 o(^tv ct .)| 2 - 


D « i (|€-...-P J --l| a -|{-...-P i | a )^-\ ( \t-...-P Nor -l\ 2 -\£<*P Ncr I 2 ) rN ° r J N « 
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0 ir;(|g--..-Pj| 2 -|g-...-Pj—W;| 2 ) * P 


n 


(|£“PAT cr . I 2 - I^PiV^I 2 ) 


Now, we freeze r 1 > r 2 > ... > , integrate out the other time variables and send e —>• 0 to 

obtain 


Jm,n* ( a ’ F a F ) 


N cr -1 


K-l r~ k [ dv [ dw f[ -± 

•/A WcP (t) +=1 l 2?r 




iV cr -l 


cr -i 

n n tv 


0(P,) - is, (If - ... - P,_, | 2 - \( - ... - P,f )/2 

R(wi) 


x 


,=o (2?r)d 0( ^ } - - • ■ ■ - Pj\ 2 - l£ - • • • - Pj - ^l 2 )/2 

1 R(£ — Po — . . . — P]y cr -l) -r-r 1 P(u>z) 


(2vr) d £!(£-•••- -PjVcr-i) - * s iV cr |£ - ■ ■ • - Pn cv - i| 2 /2 


n 


ZG^jv C7 .UB]v cr 


(2vr) d g(u; z ) + ir/|u; ; | 2 /2 




i=o 


(IAI + 1^1)! 


Here, we have changed the notation jv i->- Vi, with vq = t,VN cr+ i = 0. Next, we integrate out u>i 
except for Pi,, Pn ct , so that 

Jm,n* i' 7 ’ Fa ,£:F) 

dv f dP 

N cr (t) JlSL N crd 



tt 1 _i_ m) _ 

f = [ (2vr) rf <7(P,) - «b(|£ - • • • - Pj-il 2 - le - • • • - Pj | 2 )/2 


iVcr-l 

x Yl ( D (£ ~ ■ ■ ■ - Pj)/2) ]Ajl (D*(£ - ... - Pj)/2) lB i\(D(0)/2) lAN cr\(D*(0)/2)\ BN cr\ 

3=0 


X 


_ m-...-p Ncr - 1) 

5(C - ■ ■ • - -PjVc—l) - «hv cr |£ - ... - Pn ct 


N c , 


V2 


1 4^0 (Pn cr ) 1 2 I j 

3=0 


(Vj - Vj - +1 )l^l+l^l 

(IAI + 1^1)! 


(4.3) 
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Therefore, we have 


limE{|T e (f,0| 2 } = E E 
£—>•0 Z —' Z ' ^-0 ’ 


£—>-0 


with 


m,n =0 a:N cr >l 


11 

£—>0 


given by the RHS of (14.3|) . It is clear that n — IV cr is even, so that (—l) n fc = (— l) fc Ncr and we 
also note that 


N c , 


k-N cr = Y J {\A\ + \B i \). 


»=o 


When those crossing edges and \Aj\, \Bj\ are fixed for j = 0,..., N cr (so the RHS of (14.31) is fixed), 
the total number of possible permutations is 


(IA'I + I^I)! 


N C r 

JJ mm ' 

Now, we can sum over all permutations when N cr is fixed, denoted by (JN cr , and integrate in P^ a 
and obtain 




a N c . 


OO CXD ,, r. 

-~Uo\\i E ••• E / dv 




=0 M jVcr (i) JR(Ncr-l)d jg 0 l^'|!|^'|! 


|-4i|,|Bi|:=0 l-4jv cr li|Bjv cr l=0 ‘ 
f JV cr -l 

; ReD^-,e - ... - ij-i) 1 R eD(£ - ... - C - • • • - P Ncr - 1) 


i=i 

r iV cr -l 


x ( II (^(C- ••• -Pj)/2)\^{D*{^- ... -Pj)/2) |Bil ] (T»(0)/2) | - Ajv - l (I4*(0)/2) |BiV c.l. 
3 =o 






After the summation, we get 

z ^ *-0 ’ 


(4.4) 


Ncr 1 


a N c 


if dvf dp n 

JAN cr (t) jR(N cr -l)d 1 


,-(vj-v j+1 )R eD(£-...-Pj) 


'Ncr-l 


(4.5) 


x e -v Ncr ReD(o) | Yl ReD(P r X - ... - Pj- 1 ) ] ReD(£ - ... - P Ncr - i, £ - ■ ■ ■ - 
j'=i 
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which can also be written as 


/» I* (N C r 

E =\\4> 0 \\l / dv dP JJ < 

^ J ^ d \f=o 


,—(vj—Vj + i)ReD(£—... — Pj) 


'N cr 

x | H R eD(Pj,£ - ... - Pi- 1 ) ] 6(t - Pi - ... - PvJ. 
U =1 


(4.6) 


Thus, we have 


hmE{|T e (t,0| 2 }=ll^olll E / dv 

£ ^ 0 


I * / N C r 

/ <ip r 

JKJV cr d \j=o 


,-(Vj-Vj+ 1 ) ReD (£- -...-Pj ) 


'N c , 


(4.7) 


x | H R eD(Pj,Z - ... - P,_i) ] S(C - Pi - • • • - PvJ = Ws(t, £)• 
vi=i 


The proof of Lemma 14.21 is complete. 

Convergence of the higher order moments 

In this section, we consider convergence of the general moments 

for arbitrary M,N £ N. By Lemma 12.11 we can write 

OO 

E {Mt, 0) N } = E . . . 9m M ,s9* ni ,e ■ ■ • 9* N , E } (4.8) 

mi,...,njv=0 

with g n ,e(t,£,) = P n , e (i, 0- As f° r the variance, we only need to consider 

limJ™, (o-jP^jP), 


£—^0 




for fixed mi, ...,njv and cr such that N c {T a ) < 2. Recall that (12.2211 gives 
(<r,F a ,Z,F) | 


I I 6 

I “mi ,...,n* N ■ 


< 


irr M N 

~ 2 k dsdu / (27r)^ d e~ 0(u ' !)|,;,_? ’ r ' l/£2 5(r(;/+ rc r )P(^)d^dtc r 

1 • y<T2fc(t) ,/Ra “ fe.lvjGJV *=1 j = l 


(vi,V r )&Tc 

As before, we denote 

A* = {sjp,..., Sj,mi}, Bj = {u jt i,.. .,u jtnj }, with i = 1,..., M, j = 1,..., N. 
The pairing T a decomposes 

{A,Bj : i = 1,. • •, M,j = 1,...,N} 


26 




into the connected components. If there exists a component of size one, that is, N s (a) = 1, then, 
as in the proof of Lemma 14.11 we have a factor of 

e~ ad / 2 \4>M/e a ), 

coming from the corresponding initial condition, which implies that 

(*T Pr,£,,F) “0 

as e —> 0 since £ ^ 0. 

Thus, we only need to consider the case when 

N s (a) = N c (a) = 2. 

For any S\, S 2 £ {Ai,Bj : i = 1 = 1,...,IV} such that S\ S 2 , the following lemma 

shows that S 1 , S 2 can not be of the same type. 

Lemma 4.3. Fix a and assume N c (a ) = 2. If there exists a pair S\ , S 2 £ {A{ : i = 1,..., M} 
or S \, S 2 £ {Bj : j = 1,..., N} such that S\ -H- S 2 , then 

=°- 


Proof. Let us assume that Si = A ^, S 2 = A^ 2 - the proof for the other case is identical. Then 
we can write 


I / dsdu 

£ Jaoi-(t) 


<72 k {t) 

x \hM,hhM,i 2 \ 


n /c 2 5 { Wl + w r )^^d Wl dw r 

{vIjV^GFct 


M N 

]^[ \hM,i\ JJ 

Z=l,z/zi,Z2 j = 1 


Since A^ -H- Ai 2 and N cr (a ) = 2, after integrating in w r , we have 


h M ,h = £ ad/2 (f>o C £ a ), and h M,i 2 = e ad/2 4>o C £ a ), 


for some variable 

P = 7 ^ 0 ) 

j 

where the range of j in the summation depends on a. Now we only need to pick some p^ j and 
change this variable so that (£ — P)e a 1 — > P , which leads to 

IhM^hMMldPhd ^ e- ad ' 2 \MP)\e~ ad/2 \M^ - P)\e ad dP = \MP)M^ ~ P)\dP■ (4-9) 

Then we perform the change of variables as in the proof of Lemma 12.21 for 


M N 

n i^ini^ 

Z=l,Z^Zl,Z2 j= 1 
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and in the end obtain 


\^mi,...,n* N { a i J~ai £, F)\ 
(jM+N 


£ 2 k 


/ dsdu / 

J[ 0,i] 2fe J R kd 


n < 

(vi,Vr)eT ai 1 


-g(wdhi-^l/£ 2 R( w l) JJ (2 7 r) _d p~0( z i)hi— Wrl/e 2 


(2vr) 


%z) 


{viiVr)€.J- t,2 


2^ 

x |^ 0 (u;)^ 0 (-i - u))| n \<j>o{wi)\dwdw. 

{vi,V r )&Pa,2 


Here, as previously, zi denotes some momentum variables - we will not need their precise form, 
while (vi,v r ) £ denotes the pairings not affected by the change of variables, and (vi,v r ) £ J~a .2 
denotes the affected pairings. Finally, J - a .2 corresponds to the affected pairings when we change 
variables for 

M N 

\hM,i\ JJ \h*Nj\, 

3 = 1 

as in the proof of aforementioned Lemma 12.21 We have also changed the notation Now, 

after the temporal integration we can apply dominated convergence theorem to obtain 


4i.n^* 7 ) F O'! -F) 0: 

due to the factor — fo). □ 

By the above discussion, the nontrivial contribution of J e m ^ n *^ (a, £, F) as e —>• 0 comes 

only from the cases when M = N, the permutation a is such that 


N s {a) = N c (a) = 2, 

and all connected components contain both type-H and type-H sets. Let E(mi,..., n* N ) be the set 
of such permutations. For a £ E(m i,.. ., n* N ), we have Ai ££ Bj, i = 1,..., M, where { 1,... , M} 
is a permutation of {1 ,...,M}. We denote the set of er corresponding to a given {1 ,...,M} 
by Erj ...,n/jy). It is straightforward to check that 

M 


.^cr, £, -^) X] ''' II ^rnj.ri? ( a mi,np : £) -^)i 




where cr mi:n * denotes the permutation of Ai U Bj which keeps Ai -H- Bj. Now, we can write 
limE {* e (t,0 M (** £ (t,Z)) N } 


£—^0 


m lr ..,n N =0 <* 

oo M 

E EE E II £— >0 mi,n i ^ m *’ n i ’ ° m i-n? ’ ^ 

0 CJ " rn M ,n ^ ^ ^ 


Ml oo 


E III E E ="W*(‘.0 M 


{1,...,M}® =1 V mi,re?=0<T m4|Tl * 
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Here, the last equality comes from Lemma 14.21 


limE{|* £ (f,0| 2 }= £ X] lim J^ >r t((T rni ,r*,F^,S,F) = Ws(t,Z). 


m,i,n*=0<r n * 
1 i 

To summarize, we have shown that 


llmE 0Y } = 1 M=NM\Ws(t,£) 
£—>0 

for arbitrary M, JVsN. The proof of Theorem 11.21 is complete. 


M 


5 The fluctuation analysis 

In this section, we prove Theorems 11.31 and 11.41 
Pointwise fluctuation 

We begin with Theorem 11.31 Recall that the corrector can be written as 

OO 

U £ (t,0=e- ad/2 Y,^n,s(t^), 

n =0 

and we have previously shown that 

lim e-^M+N)/2 E {^ mit£ . . . J = E (*> ?), 

(5.1) 


a-.N s (T a )=N c (T a )=2 

when M + N = 2 K for some K G N. Let us define 

S(m i, ...,n* N ) = {a: N s {F a ) = N C (F„) = 2}. 

The constraint 

Nsi^) = N c (F a ) = 2 

forms pairings over vertices 

{Ci : 1 = 1,...,M + N} = {Ai,Bj : i = 1 ,..., M, j = 1 ,..., N}, 
or equivalently the set 

{m;,n* : i = 1,..., M,j = 1,.. .,N}. 


(5.2) 

(5.3) 


We write 


E(mi,... ,n* N ) = |JSp(mi,... ,n* N ), 


where E p (mi,...,n^) is the set of permutations corresponding to a given pairing p over (15.31) . 
Then we can write 

&'-N s (_F ?)=N c (J r ?)=2 P ff£E p (mi r ..,n^) 
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For a given p, we assume that pairs have the form (p(Z), p(Z)) with Z = 1,..., K, where 


{p(Z), p(Z') : Z = 1,..., K} = {mi,n* : i = 1,..., M, j = 1,..., N}. (5.4) 

It is straightforward to check that 

K ( 5 . 5 ) 

= E - E 

a(p(l),p(l)) <r(p(K),p(K)) 1=1 

where <r(p(Z), p(Z)) denotes the permutation of Q U Cj such that Q «->• Cj if p(Z),p(Z) corresponds 
to Ci,Cj. Now, we have 

OO 

limE {U £ (t,0 M (K(t,0) N }= E E E 

mi,...,n N =0 p C reEp(m 1 ,...,Ti^) 
oo K 

= E E E - E n>T 0 ^o,,(d^pW’P('»'^(pm.p(o)- E “f-- F ) ( 5 . 6 ) 

mi,...,n Ar =0 p CT (p(l),p(l)) a{p(K),p{K))l=l 

K / oo \ 

=EII E E_ J™ 0 ' J Ao,p(o (CT(p(z) ’ p ( z "))’- F ^(o,p( 0 )> ea ^ F ) ■ 

p 1=1 \p(0,P(0=0o-(p(0,p(0) / 

Therefore, it is clear that we only need to compute 

OO OO 

E ESo ■ 4 ,n*(T- 7 >,e Q £,i ?) and E E • 7 >’ e °£> F ) 

m,n=0 (7 m,n=0 a 


to obtain lim e _»o E{ZL e (i, £) M (ZY*(Z, £)) N }. The following lemmas combine to conclude the proof of 
Theorem 11,31 

The first lemma deals with the “complex-conjugate” moments. 

Lemma 5.1. We have 

OO 

E E lim n J Z t , n *(<r,J 7 <r,e a Z,F) = Ws,s(t,0). 

Z - ' Z ' £->0 ’ 

m,n =0 a 


Proof. Following the proof of Lemma 14.21 with £ replaced by £“£, we obtain 


OO OO ,, 

E E lim n J ^*( ( 7 ’- F -^ £ ^)^ll^ 0 ll 2 E / dv dP 


m,n =0 c 
/ N cr 

x 1 ri e 

u=° 


' N c , 


-(vj-Vj +1 )ReD(-P 0 -...-Pj) 


[] R eD(Pj, -P 0 - ... - Pj_i) ] - ... - PnJ- 

\i = 1 
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The RHS equals to Ws }S (t,0), which completes the proof. □ 

The second lemma address the “non-conjugated” moments. 

Lemma 5.2. 

OO 

E Y, lim n J kn^ :F ' T ,e a Z, F ) = W a (t,0- 

m,n =0 <t 


Proof. We use the same notation in the proof of Lemma 14.21 Recall that 
J^ n (a,F a ,e a £,F) 

= 7 T wL+i I f dsdu f TT ( 2 ' K )~ d e ~^ Wl ^ Vl ~ Vr ' [/£2 5(wi + w r )R{wi)dwidw r 
(lE) m + n J a 2 k (t) Jm. 2 kd , 


(vi,v r )&T c 


x e 


iGm (£ a £,sl m i )/e 2 e iG n (£“£,«(") ,</(") )/e 2 £ £ Pi • • • Pm ^ ^ g £ 9l ■ • • ^ 


We only need to consider a such that the number of crossing edges N cr > 1. For each crossing 
edge {rf , ), i = 1,... ,N cr , we denote the p —variable by Pj. After the integration of the delta 

functions, we obtain 


Jt, n (a,F a ,£ a ^,F) = . ^ [ dsdu [ 

m ’ nK J (*e) m+n Ja 2k (t) Jm 


dw 


fc + — AT + — 

x TT R ( w j) -8(t»j)y- fT e *(|£ a e-sR^o+-..+P-i)l 2 -h“«-Si(Po + ...+P i )| 2 )^4_ 

II ( 2vr )d II 

J=1 V ’ J=1 


^cr X 


x ii ri e 

j=0 \IgAj 
N C r 

x JJe“ |JJ 
j=i 


i(|£ Q f---Pj| 2 -|e“C---T)-w;| 2 )lL^- ] f J{\£°‘£+...+P j \ 2 -\e°‘(,+...+Pj-w l \ 2 )^d 2 


n 


V j 1 — U, 


-2(lP 7 | 2 +2P,-(P 0 +...+P,-i))A r 1 l( £a £ ~ P ° ~ ■ • • ~ PjV ^ + Pq + .. ■ + P/V e , 

cad r a ' 'rv\ a 


Compared to (14.21) . the key difference is that we get an extra factor with a large phase: 


Nc 


o -j(|P i | 2 +2P i .(P 0 +...+P i _i))r-/P 


\e 1 “* . ---^' 3 


3 = 1 


To get rid of the factor e ad , we change the variable 


(5.7) 


r n. ^ ^ — Pq — ■ ■ ■ — P N „-i + £ a Pi 


N c ., 
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Rewriting the terms in (15.71) associated with Pn ct using the new variable gives 
f, F) = , ■ \ 2k [ dsdu [ dw 

(*£) J troutt) JRkd 


^m.n ^~cr •> £ 


n 

j:Wj^P Nc , 
N cr -1 

x n > e 

3 = 1 


R(~Pq - ■■■- Pn „-1 + £ a PN cr ) - s (-P 0 -...-P Nr , r ._ 1 +£ a P N ^)^ 
{2n) d ( 2n) d 


i(| £ “?- s RPo+...+P j _i)| 2 -| £ “«-s J (Po+...+P 7 )| 2 )-W 


x e 


Rh“«-SiV c RPo + ...+PiV CT -l)| 2 -|£ a ?-SiV CT £ Q PiV c J 2 )^%^ 


N cr -1 

x n 

(TT 

.._P 3 .| 2 _| e a^..._P 3 .-^| 2 )fL_^L_ 1 r e i(| £ ^ + ... +Pj f _| £a5+ ... +Pj _^|2)h!_^i_ | 

3=0 

\l£Aj 

/ 

x n 

e i(|E“e-E“Pjv CT . 

| a -| £ “f- £ “Ply cr -«; I | a )^?- e i(h a e+£“PiV c J 2 -b“e+ e a PiV cr -^| 2 )^p- 


/Ev4iv cr 

N C r-l 

x fj e 

j=i 


-<(|P J -| 2 +2P J .(Pb+...+P i _i))4 


laB N cr 

^MZ-PnJMS + PnJ 


x e 


—z(|—Po —— -PaTct- —l+^^-PiVcr | 2 +2( — Pq — — ^Ncr — l+£ a PjV cr )’(-Po~K--+-Pv cr — l)) " 


(5.8) 

If we freeze r x ,..., r N ,P\,..., Pn ct , integrate out the other variables, and send e —>• 0, we see 
that 


lim \J^ n (a,F a ,s^,F) - H^ n (a,F a ,e^,F)\ = 0, 
£—>0 ’ ’ 


with 


F a ,e a ^ F ) = [ dv [ dP 

(-1) JA Ncr (t) jR N c d 

V cr -1 

n D(Pj, -P 0 - ... - Pj_i)/2 I D(-P 0 - • • • - -Po-...~ P Ncr - 1))/2 


I=i 

^iv cr -l 


n (^(^o+• • •+1 (^(o)/ 2 )i- 4 ^i+is Wcr i 


(5.9) 


(ujv er - UAr cr , + i) l ' 4jv ^l + |Hjv ^l jj ' ^ 


-i(\Pj I 2 +2Pj ■ (Po + ■ ■ -+Pj- i)) VJ - ■ 


(l^lVcrl + |5iV cr |)! 


- ^W cr )<M£ + -fW,,, 


i=i 


x e 


—i(|—P q —— -PA7 cr —I+^^-PaIct- | 2 +2(—Pq —— PiV cr — l+^^-PVcr )'(-P) + -”+-Pv cr — l)) " 


Here, we used the property 


D(0 = D(-d) and D(p,0 = D(-p,-£). 
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We will consider separately the cases N cr > 2 and N cr = 1. 

Multiple scattering N cr > 2. When N cr > 2, we have at least one oscillatory phase in (15.91) . 
since 


N cr ~ 1 

|| e ~ *(|fj| 2 +2Pj-(-Po+-"+-Pj-l)) 
3 = 1 


(5.10) 


x e -i(\-P 0 -...-PN c r-l+^PNcr\ 2 +^-P0---PNcr-l+^PNcrHP0 + -+PNcr-l))^T L = 


with |X| = 1 and independent of v\. For the integral in v , we have 


/» Ncr 

/ dvTUvj -v j+ i) 

J A Nr jt) f-n 


l- 4 il+l e jL _ *l- Pl l 2 72 - 


XI 


^iV er (i) j = 0 

Nc, 


p ncr pt 

<C f[ dv j\ (t-v 1 ) lAol+lBol {v 1 -v 2 ) lAll+]Bll e~ ilPll2 ^dv 1 \ 

JAjv cr -i(i) J=2 3V 2 

for some C. Applying the Riemann-Lebesgue lemma gives 

rt 


[ (t- r;i)l- 4 °l + l B °l(?;i - v 2 ) lAll+{Bll e~ ilPll2 P d Vl \ -»■ 0, 
J V2 


provided that P\ ^ 0. Thus, by the dominated convergence theorem, we obtain 

Nc 


^N cr (t) j=0 

when Pi ^ 0, which implies 


p cr 

/ dvf\(vj - v j+ i)l^l+l^l e _i|Pl|2 PX -t 0, 
J Aw_(t) Tin 


(5.11) 


(5.12) 


(5.13) 


^m,n( cr 5 Xo-, e“£, P) -» 0, as e —>• 0, 


if lV cr > 2. 


Single scattering N cr = 1. When N cr = 1, (15.91) simplifies to 


^ n (a,P CT , £ %P) = ^M f dv [ d PM^-P)M^ + P) e -^ ap I 2 

2 (- 1 ) Jo Jr 1 * 

/P(0) \ I-A>I + I®o| + |.4i| + |Bi| ^ _ v y^o| + |Bo| ^|Ai| + |Bi| 

x J WmPW 


(5.14) 


If a £ (0,1), we have a large phase factor e *l-P| 2 v/e 2 2 “ ; so f or the same reason as for N cr > 2, we 
have 


which implies 


H^ n (a,P a ,e a ^,F) ^ 0, 

d m.n( (T -d : : 7i £a iiF) v 0. 
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If a = 1, we have 


H^ n (a,P a ,E a £,F) = ^^ [ dv [ dPMi-P)^ + P)e~ i 

2 (- 1 ) Jo Ju d 

/ D(0) \ IA)|-HSo|+|.<4i|+|6i| ^ _ w y^4o|+|Bo| -yl-4i| + |Bi| 


— i\P\ 2 v 


(5.15) 


V 2 ) 


04)1 +W (|^i| + |6i|)!’ 


which is e—independent. Following the argument we used in the proof of Lemma 14.21 we have 


lin A J m,ni.^JF a ,£ a i, F ) = 11 ™ H m,n(?>Fo,e a £, F ) 

z —' £ —’ z —' >-n ’ 


cr: AT cr = l 


£—>-0 


= - D( 0, 0)e- D (°)' f dv [ dPMZ - ^)0o(£ + i>)< 

•70 7R d 

Finally, if a > 1, similarly, we have 

Y1 lim n J mA a ' ■ 7 >> e °f > F ) = 5Z lin n •?>, £% 

z —/ e—>0 z —' r— s.n ’ 


(5.16) 


—i\P\ 2 v 


a: N cr =l 


€—^0 


= - D( 0, 0)e~ D ^ t t / dPMH - P)Mt + P ). 


(5.17) 


The proof of Lemma 15.21 is complete. 

Remark 5.3. The proof shows that only single scattering contributes to the “non-conjugated” 
moments when a > 1. This is similar to the result obtained for heat equation [T] Theorem 2], 
where the single scattering constitutes the whole random corrector. For Schrodinger equation, the 
situation is different, as multiple scatterings show up in “complex-conjugated” moments as in the 
proof of Lemma 15.11 


Correlation of the fluctuations 

Here, we prove Theorem 11.41 Recall that we look at the behavior of 

W,(t,x,i) = f (5.18) 

J R d 2 2 (2vr)“ 

To prove the convergence of 

(W e (t),<p) = [ W e (t,x,€)<p*(x,£)dxd£ 

JR d 

in probability, it suffices to show the convergence of 

E{(W e (t),ip)} 

and 

E{\(W £ (t),v)\ 2 }. 
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Given that 




(27T) a y R 3d 2 2 


and 


. W6d Z Z Z Z 


x e* 71 ■* V (® 1, f 1) e~ im X2 <P(X 2 , 6) 

we first prove the following two results. 

Lemma 5.4. If a + f3 = 2 and a £ (0, 2], then as £ —>• 0, 


drjdxdf 

(27r) 2rf 


• E 

TVcr — 


e r (^ cr \ -^cr 

/ dv dP TT e -(^-^+i)^eD(-P 0 -."-P,) TT RejD (p i5 _p c 

i -^ JVcr (t) ^ R ' Vcrrf y J= o y \J=1 


-Pn - ■ ■ • - P- 




N c , 


3 =1 


X «5(--Pl - - PvJ n ^( 0 , 2 ) ll^olli + la =2 / 0 O (£ + ^ ~ ?)&(£ " “ P^P 


Lemma 5.5. //£ 1 /^ 2 , ct + /3 = 2 and a £ (0,1), then 

limER(t,6 + - £ y 1 K(*,«2 + - ^)} 

= limE{M e (t,6 + ' ^)}E{W,'((,6 + ^)M e (t,& - ^)}. 


The assumption a + f3 = 2 in Lemmas 15.41 and 15.51 matches the kinetic scaling. To see this, 
recall that 

Ue{t,0 = £ _ad/2 (V’e(L£) - ElV’eOP)}), 

and 

Mt,0 =e a ^(i/e 2 P“Oe i|£ “ ?lW - 

If we let 

^ (i, x) = cj)(t, x ) — E{c/>(f, x)}, 
then the Wigner transform written in physical domain is 



2e 


V )W*(\ 


r«+/3 


+ -ILyt-ydy, 

2e a> 


that is, we need a + (3 = 2 so that the propagation speed is of order one. Note the compensated 
phase factor from the compensation 


i£rit£ 2a +P~ 2 

^ J 
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disappears in the limit when choosing a + /3 = 2. 

Proof of Lemma 15.41 We will use the representation 

U £ (t 1 i)=e~ ad/2 Y J ^n,e{t^), 

n> 1 


so we only need to consider 


with 


E{e 


—ad cz 

17 


1 )}, 


6 = £ + 


e^rj 


e-i = e - 


£'-Ij 


2 2 

Compared to (14.21) . we need to change £ to e“£i or (the factor e“ comes from the fact that 

we are looking at the low frequency regime). Using the notations in the proof of Lemma 14.21 we 
obtain 


lim Y, ne- ad ^rn,e{t^)^ e {t^)} 
£—>•0 L ' 


m,n> 1 


/_]\n r r 

lim -—-—;— / dsdu / dw TT ' ", e 

f-vn J Lu 11 (9^r\d 


e —.0 f-—' (is) 

cr:A^ C r>l ^ 


'o- 2 fc(t) 


/i (2vr) c 


iV c , 




x n e 

3 = 1 
iVcr- / 

*n n 

i=o 

x t^o( 6---)0 o (e-i--5— 


(5.19) 


N C ; 




iPn 'TIV ■ E a+ P 2 


-j- glPj 'V r j £' 

3 = 1 


Apart from the change £ i—>■ e a £±i, the key difference between (I5.19P and (14.21) is the extra phase 
factor 

N C r 

J“J e iPj"nr-£ a+p - 2 
3 =1 

due to 7] / 0. Since a + f3 = 2, this phase factor becomes 

Ncr Ncr 

e iPj"nrJe a+f) - 2 ^ IT J p j-nr 


n 

3 = 1 


n 

3 = 1 
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and we only need to follow the proof of Lemma 14.21 to obtain 


m,n> 1 

oo « r (^ cr 

s /. .j°L/ p ( n 


-(vj-Vj + i)ReD(-P 0 -...-Pj) 


N cr =l 


■ (t) 


J ReD(Pj, -P 0 - ... - P :) 


j -U 


iVc, 


j=i 


X 5(-Pl - - P/V cr ) n W(0,2) \\M\t + la=2 / <M£ + ^ ~ p)&(£ ~ ^ ~ p)dp 


(5.20) 


The last factor comes from 


/ 

.7R' 


<M£ + £ ~y - p)4>o(£ ~~y~ p) d P . 


and the assumption of /3 = 2 — a. This finishes the proof. □ 

Proof of Lemma 15.51 The proof is similar to the case when we show the convergence of 

E{U £ (t,0 M (K(t,0) N } for M,N € N. 

The only difference is that £ is replaced by £i ± and £2 ± First, by following the proof 
of (15.61) . we have 

lim E{Z4(t,£i + ^)W*(i,£ 1 - ^)Z£(t,£ 2 + ^R(L£ 2 - ^)} 

= linrER(L£i + ^TOL£i - ^)}E{W*(t ,£ 2 + ^R(f ,£ 2 - ^)} 

+ limE{Z4(t,£i + ^)Z£(t,& + ^)}E{Z4*(t,6 - ^R(f,£ 2 - ^)} 

+ UmE{Z4(t,£i + ^R(t,£ 2 - ^)}E{Z4*(t,£ 2 + ^)Z£(t,£r - ^)} = ii + / 2 + / 3j 

e —>0 2 2 2 2 

and to complete the proof, we only need to show I 2 = / 3 = 0 . 

To study the limit of / 2 , we take, for example, 

ER(l,£i + e%/2 )«*(t,£ 2 + £^ 2 / 2 )}. 

We may follow the proof of Lemma 15.41 and obtain a phase factor 

Ncr 

J~J e iP j-(tl-&+£ l3 (ri 1 -V2)/2)r~ e “ -2 

3 = 1 

as in (15.191) . Since £1 7 ^ £ 2 , the assumption that a € ( 0 , 2 ) ensures that we have a large phase for 
multiple scattering; for single scattering, after change of variable Pi 1 —>• e“Pi, we get a factor 


c pT(£i-6+e /3 0n-^2)/2)r. e 2 “ 2 


37 






















so we have a large phase if a € (0,1). In the end, we only need to follow the proof of Lemma 15.21 
to conclude that I 2 = 0 . 

For 7; 3 , take, for example, 

E{W e (t,£i +e / V/2)Z4(t,6 -e^m/2)}. 


As in the proof of Lemma 15.21 the corresponding phase factor becomes 

Ncr 

e -i(\ p j\ 2 - p j-(£ a £i-£ a &+£ a+P {m+V 2 )/ 2 - 2 (Po+-+Pj-i))rj' /e 2 

3 =1 

as in (15.71) . The rest of discussion is the same, that is when a£ (0,1), there is always a large phase, 
which implies I 3 = 0 . □ 

Now we can discuss the limit of W £ . We use P x . to denote the Fourier transform in x,£ 
variable respectively. First, by the dominated convergence theorem, we have 

limE{(W e (t),<p>} = —i-j [ + (5.21) 

£->•0 (27T) a J R 2d £->0 2 2 

Using Lemma 15.41 we need to discuss the following two cases. 

Case 1: a + /3 = 2, a G (0,2). Using (15.201) . we integrate 77 , £ in (15.211) to obtain 


OO ft ft 

limE{(W e (t),^)} = ll^olll E / dv dP 


'N r , 


JJ e -(yj-Vj +1 )KeD{- p 0 -...-P j ) 

U=° 


N cr =l jA N cr (t) JRNard 

N cr 

R eD(Pj, -P 0 -...- Pj- 1 ) 
3 = 1 


N c , 


X s{-Pi ----- PNjp&n- E p i y p °) 

3= 1 


Ws j 3 (t,x, 0 )(p*(x,€)dxd£, 


with 


'N c , 


p p 1 2.1 cr 

w SjS (t,x,o= WfaWl V / dv dP \] 


,-(yj- v j+i)ReD(£-Pb--.-Pj) 


3=0 


Ncr 

x | ReD(Pj,£ - P 0 
U =1 


Nc, 


■ . - Pj—i) ] 6(£ - Pi - ... - PncMx - + E PjVj). 

3 = 1 


Clearly, we have 

w s>8 ( t,*,o = w^(t,s,o - n<A)ii2‘J(e)<y(*)e- Ra2J(0)t , 

which consists of the scattering component of the transport equation (11.281) with the initial condition 


Ws( 0 ,x ,0 = ll^oll^^)^^)- 
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Case 2: a = 2, f3 = 0. By a similar discussion, we have 


lim E{(W £ (f), </?)} = / W StS (t,x,0)<p*(x,£)dxd£; 

0 ./iu 2 d 


(5.22) 


with 

Ws,s{t,x,Z) = Ws(t,x,0 - (2TT) d 5(O\Mx)\ 2 e- ReD{0)t , 
and Ws(t,x, £) solving (|1.28[) with initial condition Ws(0 ,x,g) = (27r) rf <5(^)|0o| 2 • 
By Lemma 15.51 if we further assume a £ (0,1), we have 

limE{l(W £ (t),p)l 2 } = llimE{(W £ (t),<p)}l 2 , 

E-tO e-tO 

which implies (W £ (t),ip) converges in probability. 


(5.23) 


A Moments of product of Gaussians 

The following result is standard, we present a proof for the sake of convenience. We assume that 

{Nij : i = 1,. ,.,m, j = 1,..., M t } 

are zero-mean real (complex) Gaussian random variables, and write 

m Mi 

E (nn^}=E n ^ 

4=1 J =1 T ((ij'),(iJ))G^ 

where Ylr extends over all pairings formed over vertices {(i. j) : i = 1,... ,m,j = 1,..., Mj}. We 
set 

A = {( i,j) : j = 1 , • ■ •, Mi}. 

For a given pairing F and i ^ i, we say that A % is connected to A, and denote this by A. L A, if 
there exist j, j such that (( i,j ), (i,j)) € T. In this way, the set {A : i = 1,... ,m} is decomposed 
into connected components, and we denote the size of the smallest component by N S (F). 

Lemma A.l. For each i = 1,..., m, let Xj = fljii Aj> then we have 

m 

E{[[(X i -E{X i })}= E I! ( A - 2 ) 

*=i f -. n s ( t )>2 

Proof. We write 

mm m 

- E {XJ)} = E {A - E {^})} - E{A!}E{[](W - E{W})}, (A.3) 

i =1 »=2 i =2 
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and note that every term in the expansion of 


nix, - E{Xi}) 


i=2 


is a product of zero-mean Gaussians (with possible multiplicative constant), so when taking expec¬ 
tation, we follow the rule of computing joint moments of zero-mean Gaussians. For any pairing 
such that A\ is not connected to any A t , i ^ 1 , we have a cancellation from the corresponding term 
in 


EpGjEfQpQ-ElxJ)}. 


i=2 

Thus, we can write 

m m 

E{J](*i - E {^})} = H(Xi - E{XJ)}, (A.4) 

i= 1 i =2 

where Ei stands for the expectation with the summation over those T such that A\ -H- Ai for some 
i / 1. Following a similar procedure for X 2 — EjA^}, we have 


Ein^ - 


(A.5) 


i =1 


i=3 


with Ei ; 2 stands for the expectation with the summation over those T such that A\ Ai for some 
j / 1 and A 2 <->■ Ai for some i 7 ^ 2. In the end, we obtain 


E{f](^- E {^})} = E i„.., m Alx,h 


(A. 6 ) 


i =1 


i =1 


where we only take the expectation with the summation over those J- such that for alii = 1 ,,m, 
Ai Aj with some j 7 ^ i, and these are exactly the pairings with iV^X) > 2 . □ 
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